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ABSTRACT 

Recently, a soft black-body component was observed in the early X-ray afterglow of 
GRB 060218, which was interpreted as shock breakout from the thick wind of the pro- 
genitor Wolf-Rayet (WR) star of the underlying Type Ic SN 2006aj . In this paper we 
present a simple model for computing the characteristic quantities (including energy, 
temperature, and time-duration) for the transient event from the shock breakout in 
Type Ibc supernovae produced by the core-collapse of WR stars surrounded by dense 
winds. In contrast to the case of a star without a strong wind, the shock breakout oc- 
curs in the wind region rather than inside the star, caused by the large optical depth 
in the wind. We find that, for the case of a WR star with a dense wind, the total 
energy of the radiation generated by the supernova shock breakout is larger than that 
in the case of the same star without a wind by a factor > 10. The temperature can 
be either hotter or cooler, depending on the wind parameters. The time-duration is 
larger caused by the increase in the effective radius of the star due to the presence of 
a thick wind. Then, we apply the model to GRB 060218/SN 2006aj. We show that, 
to explain both the temperature and the total energy of the black-body component 
observed in GRB 060218 by the shock breakout, the progenitor WR star has to have 
an unrealistically large core radius (the radius at optical depth of 20), larger than 
lOOi?©. In spite of this disappointing result, our model is expected to have important 
applications to the observations on Type Ibc supernovae in which the detection of 
shock breakout will provide important clues to the progenitors of SNe Ibc. 
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1 INTRODUCTION 



direct observations of supernova features in the spectra of 
GRB afterglows, or from indirect observations of rebright- 
ening and/or flattening (called "red bumps") in GRB af- 
terglows which are interpreted as the emergence of the un- 
derlying supernova lightcurves (Delia Valle 2006; Woosley 
& Bloom 2006; Woosley & Heger 2006b, and references 
therein). The discovery of the connection between GRBs and 
supernovae has been one of the most exciting developments 
in the fields of GRBs and supernovae in the past decade. 



Since the first detection of the afterglows (Costa et al. 1997; 
van Paradijs et al. 1997; Frail et al. 1997) and the host 
galaxies (Bloom et al. 1998, 1999; Fruchter et al. 1999b) of 
gamma-ray bursts (GRBs), by now it has well been estab- 
lished that long-duration GRBs are cosmological events oc- 
curring in star-forming galaxies (Paczynski 1998a; Fruchter 
et al. 1999a; Berger, Kulkarni & Frail 2001; Frail et al. 
2002; Christensen, Hjorth & Gorosabel 2004; Sollerman et 
al. 2005; Fruchter et al. 2006, and references therein), and 
are most likely produced by the core-collapse of massive 
stars (Woosley, Heger & Weaver 2002; Piran 2004; Zhang 
& Meszaros 2004; Woosley & Heger 2006a, and references 
therein). This scenario has received strong support from 
the cumulative evidence that some, if not all, long-duration 
GRBs are associated with supernovae (SNe), either from 



Interestingly, all the supernovae that have been spec- 
troscopically confirmed to be associated with GRBs, in- 
cluding SN 1998bw/GRB 980425 (Galama et al. 1998), 
SN 2003dh/GRB 030329 (Hjorth et al. 2003b; Stanek et 
al. 2003), SN 20031w/GRB 031203 (Malesani et al. 2004; 
Sazonov, Lutovinov & Sunyaev 2004), and the most recent 
one, SN 2006aj/GRB 060218 (Masetti et al. 2006; Modjaz et 
al. 2006; Campana et al. 2006; Sollerman et al. 2006; Pian et 
al. 2006; Mirabal et al. 2006; Cobb et al. 2006), are Type Ic 
having no detectable hydrogen and helium lines. However, 
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the supernovae that are associated with GRBs also remark- 
ably differ from ordinary Type Ibc supernovae: they have ex- 
tremely smooth and featureless spectra indicating very large 
expansion velocity, are much more energetic (i.e., involving 
much larger explosion energy), and eject significantly larger 
amount of nickels (Hamuy 2004; Delia Valle 2006; Woosley 
& Heger 2006b), except SN 2006aj/GRB 0602f8 which is 
somewhat closer to normal SNe Ibc (see below; Mazzali et al. 
2006). For these reasons, they are often called "hypernovae" 
to be distinguished from normal supernovae (Iwamoto 1998; 
Paczyhski 1998a,b). A correlation between the peak spec- 
tral energy of GRBs and the peak bolometric luminosity of 
the underlying supernovae are found by Li (2006), based on 
the multi-wavelength observations on the above four pairs 
of GRBs-SNe. 

The discovery of GRB-SN connection has provided us 
with important clues to the progenitors of GRBs, since it 
is broadly believed that Type Ibc supernovae are produced 
by the core-collapse of Wolf-Rayet (WR) stars who have 
lost their hydrogen (possibly also helium) envelopes due to 
strong stellar winds or interaction with companions (Smartt 
et al. 2002; Woosley et al. 2002; Filippenko 2004; Woosley 
& Heger 2006a, and references therein). In fact, for several 
GRBs, observations with high quality optical spectra have 
identified the presence of highly ionized lines with high rel- 
ative velocities most likely coming from shells or clumps of 
material from WR stars, supporting WR stars as the GRB 
progenitors (Mirabal et al. 2003; Schaefer et al. 2003; Klosc 
et al. 2004; Chen, Prochaska & Bloom 2006, see, however, 
Hammer et al. 2006). 

A systematic study on the GRB afterglows carried out 
by Zeh, Klose & Hartmann (2004) suggested that all long- 
duration GRBs are associated with supernovae. However, it 
appears that only a small fraction of Type Ic supernovae are 
able to produce GRBs, since the rate of GRBs and hyper- 
novae are several orders of magnitude lower than the rate 
of core-collapse supernovae (Podsiadlowski et al. 2004). Al- 
though both long-duration GRBs and core-collapse super- 
novae are found in star-forming galaxies, their location in 
the hosts and the morphology and luminosities of their host 
galaxies are significantly different as most clearly revealed 
by the recent study of Fruchter et al. (2006) with Hubble 
Space Telescope (HST) imaging. The core-collapse super- 
novae trace the blue-light of their hosts that are approxi- 
mately equally divided between spiral and irregular galaxies, 
while long GRBs are far more concentrated on the brightest 
regions of faint and irregular galaxies. Fruchter et al. (2006) 
argued that their results may be best understood if GRBs 
are formed from the collapse of extremely massive and low- 
met allicity stars. 

The preference of long-duration GRBs to low- 
metallicity galaxies (Fynbo et al. 2003; Hjorth et al. 2003a; 
Le Floc'h et al. 2003; Sollerman et al. 2005; Fruchter et al. 
2006) has been strengthened by the recent paper of Stanek 
et al. (2006), in which a strong anti-correlation between the 
isotropic energy of five nearby SN-connected GRBs and the 
oxygen abundance in their host galaxies was found, which 
was used to argue that the life in the Milky Way is protected 
away from GRBs by metals. Stanek et al. (2006) have sug- 
gested that long-duration GRBs do not trace star formation, 
but trace the metallicity. 

The discovery of GRB 060218 and its association with 



SN 2006aj by Swift has shed more light on the GRB-SN 
connection as well as on the nature of GRBs. GRB 060218 
has a cosmological redshift z = 0.0335 corresponding to a 
luminosity distance of 147Mpc (f2 m = 0.3, SIa = 0.7, and 
Ho = 70 km s~ 1 Mpc _1 ), which makes it the second near- 
est GRB among those having determined redshifts (about 
four times the distance of GRB 980425 at z = 0.0085; Cam- 
pana et al. 2006; Pian et al. 2006; Sollerman et al. 2006). 
GRB 060218 is very unusual in several aspects. It has an 
extremely long duration, about 2, 100 s. Its spectrum is 
very soft, with a photon index 2.5 ± 0.1 and peak energy 
Sp 0a k = 4. 91q 3 keV in the GRB frame. The isotropic equiv- 
alent energy is Ei BO = (6.2 ±0.3) x 10 49 ergs extrapolated to 
the 1-10, 000 keV in the rest frame energy band (Campana 
et al. 2006), which is at least 100 times fainter than nor- 
mal cosmological GRBs but among a population of under- 
energetic GRBs (Sazonov et al. 2004; Liang, Zhang & Dai 
2006). 

Although the supernova associated with GRB 060218, 
i.e. SN 2006aj, is broadly similar to those previously dis- 
covered GRB-connected supernovae, it also shows some re- 
markable unusual features (Pian et al. 2006; Sollerman et al. 
2006; Mazzali et al. 2006). Among the four GRB-connected 
supernovae mentioned above, SN 2006aj is the faintest one, 
although still brighter than normal Type Ibc supernovae. 
Its lightcurve rises more rapidly, and its expansion velocity 
indicated by the spectrum is intermediate between that of 
other GRB-connected supernovae and that of normal SNe 
Ibc. Modeling of the spectra and the lightcurve of SN 2006aj 
reveals that SN 2006aj is much less energetic compared to 
other GRB-connected supernovae: it had an explosion en- 
ergy E- ln ~ 2 x 10 51 ergs, ejected a mass M e j « 2Mq, com- 
pared to E in ~ 3-6 x 10 52 ergs, and M cj ~ 1OM of the 
others (Mazzali et al. 2006). This suggests that SN 2006aj 
is closer to normal Type Ibc supernovae than to the other 
GRB-connected supernovae, and there does not exist a clear 
gap between hypernovae and normal Type Ibc supernovae 
(Li 2006). 

The X-ray afterglow observation by the X-Ray Tele- 
scope (XRT) on board Swift on GRB 060218 started 159 s 
after the burst trigger. A very interesting feature in the early 
X-ray afterglow is that it contains a soft black-body compo- 
nent which has a temperature about 0.17 keV and comprises 
about 20% of the total X-ray flux in the 0.3-10 keV range, 
lasting from 159 s up to ~ 10, 000 s. The black-body compo- 
nent was not detected in later XRT observations (Campana 
et al. 2006). The total energy contained in the black-body 
component, as estimated by Campana (private communica- 
tion), is « 10 49 ergs. Campana et al. (2006) interpreted it 
as supernova shock breakout from a dense wind surrounding 
the progenitor WR star of the supernova. 

Butler (2006) conducted an analysis on the early X- 
ray afterglows of a sample (> 70) of GRBs observed by the 
XRT / Swtft. He found that although most of the afterglow 
spectra can be fitted with a pure power law with extinc- 
tion, a small fraction of them show appreciable soft thermal 
components at 5-10% level. His reanalysis on GRB 060218 
showed that the black-body component contains energy as 
much as 2.3 x 10 50 ergs and has a duration ~ 300 s. Accord- 
ing to Butler's analysis, the soft black-body component even 
dominates the flux after ~ 1, 000 s from the burst trigger. 

Flashes from shock breakout in core-collapsed super- 
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novae were first predicted by Colgate (1968) almost forty 
years ago, originally proposed for GRBs that had not been 
discovered yet. However, they have not been unambigu- 
ously detected in supernova observations yet (Calzavara & 
Matzner 2004). This is mainly due to the transient nature 
of the event. It is generally expected that the flash from 
shock breakout precedes the supernova, is much brighter 
and harder than the supernova radiation but has a very 
short time-duration. 

According to the general theory of core-collapsed su- 
pernova explosion, the liberation of explosive energy in the 
interior of a progenitor star generates a shock wave. The 
shock wave propagates outward. However, the external ap- 
pearance of the star remains unaltered, until the shock wave 
reaches a point (the shock breakout point) near the stellar 
surface where the diffusion velocity of photons begins to ex- 
ceed the shock velocity. The postshock radiation can then 
leak out in a burst of ionizing radiation, producing a bril- 
liant flash in the UV/X-ray band (Klein & Chevalier 1978; 
Chevalier & Klein 1979; Imshennik & Nadezhin 1989, for a 
comprehensive review see Matzner & McKee 1999). 

For the famous Type II SN 1987A, theoretical calcula- 
tions have shown that the shock emergence from the sur- 
face of the progenitor (Sk 1, a blue supergiant) would have 
produced a radiation of ~ 10 47 ergs in the extreme UV to 
soft X-ray band, lasting 1-3 minutes (Imshennik & Nadezhin 
1988, 1989; Ensman & Burrows 1992; Blinnikov et al. 2000). 
In fact, in the observed bolometric lightcurve of SN 1987A, 
there was a fast initial decline phase which could be the tail 
of the lightcurve produced by the shock breakout (Imshen- 
nik & Nadezhin 1989). If the shock breakout interpretation 
of the soft black-body component in GRB 060218 is con- 
firmed, it would have important impact on the theories of 
both GRBs and supernovae. The case of GRB 060218/SN 
2006aj would also be the first unambiguous detection of a 
shock breakout event from supernovae. 

Although the propagation of a strong shock in a super- 
nova and the appearance of shock emergence (shock break- 
out) have been intensively studied both analytically and nu- 
merically (Klein & Chevalier 1978; Imshennik & Nadezhin 
1988, 1989; Ensman & Burrows 1992; Blinnikov et al. 1998, 
2000, 2002; Matzner & McKee 1999; Tan, Matzner & McKee 
2001), in the situation of supernovae produced from stars 
with dense stellar winds they have not been fully explored 
yet. If the stellar wind of the progenitor is very optically 
thick — which is indeed the case for Type Ibc supernovae 
whose progenitors are believed to be WR stars — the shock 
breakout will occur in the wind region after the shock passes 
through the surface of the star, instead of in the region in- 
side the star. Since a stellar wind has a mass density profile 
very different from that of a star, the model that has been 
developed for the shock emergence in supernovae with pro- 
genitors without stellar winds cannot be directly applied to 
the case of progenitors with dense stellar winds. 

In this paper, we present a simple model for semi- 
analytically computing the propagation of a strong shock 
in a dense stellar wind, and estimating the characteristic 
quantities for the transient event from the shock breakout 
in SNe Ibc. The model is obtained by an extension of the 
existing model for the shock propagation and breakout in 
supernovae produced by the core-collapse of stars without 
dense stellar winds. Then, we apply the model to SN 2006aj 



and examine if the soft black-body component in the early 
X-ray afterglow of GRB 060218 can be interpreted by the 
supernova shock breakout. 

The paper is organized as follows. In Sec. 2, we describe 
a simple but general model for the mass density and veloc- 
ity profile for the wind around a WR star, and calculate the 
optical depth in the wind. In Sec. 3, we model the propaga- 
tion of a supernova shock wave in a stellar wind, taking into 
account the relativistic effects. In Sec. 4, we analyze the evo- 
lution of the shock front, and the radiation energy contained 
in it. In Sec. 5, we present a procedure for calculating the 
quantities characterizing the transient event arising from the 
shock breakout, including the released energy, the temper- 
ature, the time-duration, and the momentum of the shock 
front at the time of shock breakout. In Sec. 6, we present 
our numerical results. In Sec. 7, we apply our model to GRB 
060218/SN 2006aj. In Sec. 8, we summarize our results and 
draw our conclusions. 

Appendix A is devoted to the formulae for computing 
the optical depth of a wind in the framework of the stan- 
dard stellar wind model. Appendix B lists the formulae for 
computing the characteristic quantities for supernova shock 
breakout from a star without winds, in the trans-relativistic 
regime. Appendix C presents a correlation in the WR star 
parameters. 



2 MASS DENSITY PROFILE OF THE WIND 
OF A WOLF-RAYET STAR AND THE 
OPTICAL DEPTH IN THE WIND 

Wolf-Rayet stars are very luminous, hot, and massive stars 
that are nearly at the end of their stellar lives. Based on their 
spectra, WR stars are often classified as WN stars (nitrogen 
dominant) and WC stars (carbon dominant). WR stars are 
characterized by extremely dense stellar winds, losing mass 
at a rate of 10 _6 -10 _4 Mq yr _1 with a wind terminal veloc- 
ity of 700-3, 000 km s -1 . The mass lost from WR stars by 
stellar winds is so enormous that most (if not all) of hydro- 
gens of WR stars have been lost. This is a main reason for 
the general belief that WR stars are the progenitors of Type 
Ibc supernovae. 

Some basic relations in the physical parameters of WR 
stars can be found in Langer (1989), Schaerer & Maeder 
(1992), and Nugis & Lamers (2000). 

The wind of a WR star is usually extremely dense. This 
is characterized by the fact that, for the majority of WR 
stars, the ratio of the momentum of the wind (Mvoo , where 
M is the mass-loss rate, Voo is the terminal velocity of the 
wind) to the momentum of radiation (L/c, where L is the 
luminosity of the star, c is the speed of light) is much larger 
than unity, indicating that on average each photon leaving 
the star must be scattered several times and the wind must 
be optically thick. As a result, the photospheric radius (R p h, 
the radius where the optical depth t w = 2/3) often differs 
from the core radius of the star (R*, the radius where t w = 
20 by definition) by a factor > 2. 

In Fig. 1, we plot the photospheric radius against the 
core radius for 86 Galactic WC-type and WN-type stars 
(Hamann, Koesterke & Wessolowski 1995; Koesterke & 
Hamann 1995) and 6 LMC WC-type stars (Grafener et 
al. 1998), determined with the "standard model" of stel- 
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Figure 1. Photosphcric radius versus stellar core radius, for 86 
Galactic WRs (filled circles for WC-type, open circles for WN- 
type) and 6 LMC WRs (triangles, WC-type only). The dashed 
lines show the relation of R p ^ = J?*, R p ^ = 2R+, and R p ^ = 10R* 
(upward). 



lar winds. For many WRs, especially those of WC-type, we 
have R p h > 2R+. 

The mass density of a steady and spherically symmetric 
wind is related to the mass-loss rate and the wind velocity 
by 



p(r) = 



M 



4nr 2 v r (r) 



(1) 



where r is the radius from the center of the star, and v r is 
the radial velocity of the wind. We model the velocity of the 
wind by (Schaerer 1996; Ignace, Oskinova & Foullon 2000; 
Nugis & Lamers 2002) 



V r (r) =Voo(l- ' 



(2) 



where a < 1 and b > 1 are free parameters. The presence of 
q in equation (2) is to ensure that the mass density of the 
wind is regular at the stellar radius r — R+. 

In the "standard model" of stellar winds the value of b 
is assumed to be unity, as in the case of O-stars. However, 
it has been argued that for WR stars b can be significantly 
larger (Robert 1994; Schmutz 1997; Lepine & Moffat 1999). 
According to the calculations of Nugis & Lamers (2002), b 
is typically in the range of 4-6. 

The value of a can be determined by the radial velocity 
of the wind at r = i?*. If we define e = v+fvoo, where v+ = 
v r (R±), then 



.1/6 



(3) 



Typically, has the order of the sound speed at 7i*, and 
e ~ 0.01 (Schaerer 1996). 

In the outer wind region, where r 3> i?* and v r ~ «oo, 
the wind density p ~ r -2 . In the region close to the stellar 




Figure 2. The log-slope of the wind density, s = dlnp/dlnr = 
—2 — b(r/aR+ — l) -1 , as a function of the radius r. As r — > 
co, s approaches —2 (the upper dashed line). For small r, s is 
significantly smaller than —2. A shock wave accelerates only in 
the region of s < —3 (below the lower dashed line; see Sec. 3). 
Left to right: b = 1-10 with Ab = 1. 



surface (i.e., r ~ R*), the wind density has a much steeper 
log-slope (Fig. 2). As will be seen in Sec. 3, it is the very 
steep mass density profile near the surface of the star that 
makes it possible for a shock wave to accelerate in the wind 
region. We will also see in Sec. 6 that shock breakout takes 
place at a radius not far from the surface of r = i?*. So 
we adopt equation (2) for the wind velocity profile since its 
asymptotic form v r = t>oo (and p oc r~ 2 ) is not accurate for 
describing the wind velocity (and hence the mass density) 
near r = _R*. 

The opacity n w in the WR wind region is complex 
and generally a function of radius (Nugis & Lamers 2002; 
Grafener & Hamman 2005). However, compared to the 
mass density, the opacity changes very slowly with ra- 
dius. For example, at the sonic point in the wind, we have 
dlnK w /dlnr ~ 0.001-0.03 (Nugis & Lamers 2002), while 
|dlnp/dlnr| > 2 always. Hence, to calculate the optical 
depth in the wind, we can approximate k w by a constant 
although its value is uncertain at some level. Then, the op- 
tical depth in the wind is 



/oo 



pdr - 



A 



(b-l)aR* 



1-° 

y 



(4) 



where y = r / R* and A = k w M /(AitVoo)- 

As commonly adopted in the literature, we define the 
stellar core radius i?* of the WR star to be the radius where 
t w — 20. Then, we can rewrite the optical depth as 



Tw = TO 



1-° 

y 



(5) 
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Figure 3. The solution for a and j/ pn = i? pn /-R*, for e = 10 5 — 
1CT 1 . Top to bottom: b = 1 - 10 with Ab = 1. The filled circles 
on each curve label the values of e = 10~ 5 , 10~ 4 , 10~ 3 , 10~ 2 , 
and 10" 1 from left to right. 



where 



TO = 



20 



1 



(6) 



(b-l)aR* (l-a) 1 - 6 

By definition, the boundary of the photosphere is at the 
photospheric radius where t w — 2/3. Then, we can solve for 
yph = Rph/R* from equation (5) 

i/(i-6)-i - 1 



(7) 



For a given 6, y p h is a decreasing function of a. As 
a — ► 1, we have to — > and j/ p h — > 1. As a — > 0, we have to « 
20/ [(6 — 1)q] and t/ p h — > 30. Thus, in general we must have 
1 < J/ph < 30. By equation (3), a is a decreasing function of 
e. 

The above results for the optical depth are valid only for 
b > 1. The corresponding formulae for 6 = 1 (the "standard 
model") are given in Appendix A. 

In Fig. 3, we plot a versus y p h, solved from equations (3) 
and (7) (or eq. A3 when 6 = 1) for a set of discrete values 
of b (from 1 to 10) and continuous values of e (from 10~ 5 
to 10 _1 ). The value of j/ph sensitively depends on e. For the 
same value of b, y pb drops quickly as e decreases. When e 
is fixed, y p h decreases if b increases from 6 = 1, very fast 
for small values of e. However, for b > 3, the effect of the 
variation in 6 on the value of j/ p h is not dramatic. 

The opacity k w , and the corresponding optical depth 
t w , are for the optical photons in the wind and are hence 
valid for calculating the mass density profile of the wind be- 
fore a supernova shock passes through it. As will be seen in 
the following sections, the radiation generated by the super- 
nova shock wave in the wind of a WR star is in the X-ray 
band and we also need consider the opacity and the optical 
depth to the X-ray photons for computing the thickness of 



the shock front and the emergence of the shock wave (Sees. 4 
and 5). 

The X-ray opacity of a gas strongly depends on the ion- 
ization state of the gas (Krolik & Kallman 1984). The radi- 
ation generated by a supernova shock wave in the wind of a 
WR star has a luminosity Lx > 10 45 ergs s _1 (Sees. 5 and 6), 
which contains enough photons to fully ionize the surround- 
ing gas. This fact can be seen from the ionization parameter, 
defined as the ratio of the photon number density to the par- 
ticle number density, S = Lx/ (47rcr 2 nH£ p h) , where e p h is 
the energy of photons. Using equation (1) (where v r « Voo) 
and uh = p/f^Hinu, where mn is the mass of proton, and 
fiu ~ 2 is the mean molecular weight per proton, we get 

oc 



Me ph c 
4.4 x 10 6 ft' 1 1 



\ 10 45 ergss~ 



(i 



£ph 



keV 



where 



fi = 



M 



5 x 10- 5 M©yr- 



2, 000 km s" 



(8) 



(9) 



Hence, for typical parameters we have E > 10 , which 
means that a tiny fraction of the radiation would be enough 
to fully ionize the gas in the wind. Then, the absorption 
opacity is negligible (Krolik & Kallman 1984). The opac- 
ity in the wind to the X-ray photons is then dominated by 
the electron scattering opacity, k cb — 0.2 cm 2 g _1 , which 
is insensitive to the photon energy if the photon energy is 
much smaller than the electron mass energy (Akhiezer & 
Berestetskii 1965). 

The X-ray optical depth in the wind is 



Tx 



f 

foes I 
J r 



pdr - 



(10) 



The X-ray photospheric radius, defined by rx = 2/3, is then 



J/ph,x = a 



+ 3(,To) 



1/(1-6) 



(11) 



3 PROPAGATION OF A STRONG SHOCK 
WAVE IN THE WIND 

The propagation of a strong shock wave in a gas is deter- 
mined by two competing processes: the collection of mass 
from the ambient gas makes the shock wave decelerate, and 
the steep downward gradient of the gas mass density makes 
the shock wave accelerate. Based on previous self-similar an- 
alytical solutions and numerical works, Matzner & McKee 
(1999) have proposed a continuous and simple form for the 
shock velocity that accommodates both spherical decelera- 
tion and planar acceleration 



/Sin 
V 3 OC 

V m 



1/2 



pi- 
rn 



-01 



(12) 



where (3i w 0.2. In the above equation, E ln is the explosion 
kinetic energy, m(r) = M(r) — M rcm , M rcm is the mass of 
the material that will become the supernova remnant, and 
M(r) is the mass of the material contained in radius r. 

After the shock has collected an enough amount of mass 
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so that m(r) does not change significantly any more, we have 
v a oc (pr 3 ) 1 , the behavior of the shock is purely deter- 
mined by the profile of the mass density in the region that 
the shock is plowing into. Then, for a spherically symmetric 
gas, the shock wave accelerates when d (pr 3 ) /dr < 0, and 
decelerates when d (pf 3 ) /dr > 0. 

To generalize the formalism to the case of a relativis- 
tic shock wave, Gnatyk (1985) has suggested to replace the 
shock velocity v s on the left-hand side of equation (12) 
by r s /3 s , where (3 S = v a /c, c is the speed of light, and 

r s = (l — 0f) 1 '" 2 is the Lorentz factor. The equation so 
obtained quite accurately describes both the limits of non- 
relativistic {01 <C 1, i.e., F s « 1) and ultra-relativistic 
(01 « 1, i.e., V a > 1) shocks. However, Tan et al. (2001) 
have shown that it is not accurate in the trans-relativistic 
regime (01 close to 1 but Y a not large enough). Tan et al. 
(2001) have suggested the following formula for both non- 
relativistic and trans-relativistic shocks 

2\0.12 



Ts0 = p(l+p 2 )' 



p = 0.736 



Vmc 2 



m 



(13) 



With numerical simulation, Tan et al. (2001) have ver- 
ified equation (13) for trans-relativistic and accelerating 
shocks with T a /3 a up to a value ~ 10. However, the lim- 
ited numerical resolution in their code has not allowed them 
to follow the acceleration of a non-relativistic shock into the 
ultra-relativistic regime (Tan et al. 2001). 

Although equation (13) has never been tested on rela- 
tivistic and decelerating shock waves, in the non-relativistic 
limit it returns to the formula of Matzner & McKee, i.e. 
equation (12), which applies to both accelerating and decel- 
erating shocks. Hence, we assume that equation (13) applies 
to both accelerating and decelerating relativistic shocks. 

Because of the compactness of WR stars, the problem 
that we are solving here is right in the trans-relativistic 
regime (as will be confirmed latter in this paper). Thus we 
will use equation (13) to calculate the momentum of a shock 
wave propagating in a wind of a WR star. In addition, since 
the wind contains a negligible amount of mass, at the radius 
where shock breakout takes place (either inside the star but 
close to its surface, or in the wind region) , we have m « M ei , 
where M e j is the ejected mass. 

Although the equation for the shock movement that 
we use in this paper is the same as that used by Matzner 
& McKee (1999) and Tan et al. (2001), the mass density 
profile in the wind of a star is very different from that in 
the interior of a star. In the outer layer of a star the mass 
density drops quickly as the radius increases by a very small 
amount, as described by equation (Bl), Hence, as the shock 
wave approaches the surface of the star, it always accelerates 



according to v s oc p Pl since m ■ 



M c j and r : 



The situation is very different in a stellar wind. A shock 
wave propagating in a wind with a density given by equa- 
tions (1) and (2) accelerates in the region near the stellar 
surface r = R*, but decelerates at large radius since p oc r~' 2 
and d(pr 3 )/dr > for r 3> R* (Fig. 2). The transition from 
acceleration to deceleration occurs at a radius determined by 
d(pr 3 )/dr — 0, where the shock velocity reaches the maxi- 



mum. The transition radius is found to be 

R a = (l + 6)ai?* . (14) 

After passing the transition radius, the shock wave 
starts decelerating. The maximum T S /3 S is then obtained by 
submitting r — R a into equation (13) 



= Pmax(l+Pm 

= 1.181 [«/(&)] 



1/2 



(15) 



where f(b) = (1 + b) (1 + l/b)\ and the mass function 'J' is 
defined by 



MR, 



1.654 x 1O" 9 M p 



SR7. 



(16) 



where p is defined by equation (9). 

The function Vl/ is an estimate of the mass contained in 
the photosphere region of the wind. A correlation between 
*!/ and -Rph is presented in Appendix C. 

Submitting fiducial numbers in, we get 

.... n -0.187 / „ \ 0.5 

<xf(b) 



p. 



= 1.137p" 



\ WM Q J 



/(5) 

0.313 



3R Q J 



E in 
10 52 ergs 

0.187 



(17) 



Thus, typically, the shock wave is trans-relativistic. 



4 ENERGY OF THE RADIATION 

CONTAINED IN THE SHOCK FRONT 

The gas pressure behind a relativistic shock front, measured 
in the frame of the shocked gas, is (Blandford & McKee 
1976) 



P2 = (72 — 1)(7272 + l)pc , 



(18) 



where 72 is the polytropic index of the shocked gas, 72 is 
the Lorentz factor of the shocked gas, and p is the mass 
density of the unshocked gas. The Lorentz factor 72 is re- 
lated to the Lorentz factor of the shock front r = T s by the 
equation (5) of Blandford & McKee (1976). Since WR winds 
are radiation-dominated, we have 72 = 4/3. Then, we can 
approximate P2 by 



P2 w F p (r s v s ) pT 2 s v 2 s , 
where F p (T s v s ) ~ 1 is defined by 
2 A 



3 21 (1 + 0.4252 x 2 ) - 4144 ' 



(19) 



(20) 



which has the correct asymptotic values as r s v s 00 (ultra- 
relativistic limit) and r s v s — > (non-relativistic limit), and 
has a fractional error < 0.3% in the trans-relativistic regime. 

Denoting the temperature of the radiation behind the 
shock front by T2, then the pressure of the radiation mea- 
sured in the frame of the shocked gas is 



iaT 2 4 ^ p 2 w F p {T a v s ) pV 2 a v 2 s , 

where a is the radiation density constant. 



(21) 
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A strong shock has a very narrow front. In the non- 
relativistic limit, the geometric thickness of the shock front 
is (Imshcnnik & Nadezhin 1988, 1989) 

Ar s « , (22) 

where n CB is the electron scattering opacity (see Sec. 2). That 
is, the thickness of the shock front is equal to the mean free 
path of photons multiplied by the optical depth of the shock 

t. = - . (23) 

V s 

For an ultra-relativistic blast wave, the total energy 
stored in the shock wave is proportional to r 2 r 3 and hence 
the thickness of the shell of shocked particles is ~ r /r 2 
(Blandford & McKee 1976). That is, in the ultra-relativistic 
limit, the thickness of the shock front measured in the rest 
frame is oc T~ 2 . Hence, using the optical depth of the shock 
and that of the wind, we can estimate the geometric thick- 
ness of a relativistic shock front in the rest frame by 



(24) 



where rx is the X-ray optical depth (eq. 10) , and 
(9hir x 



TX 
KesPV 



a(b- 



<91nr 



1--- 1 



V 



(25) 



Equation (24) returns to equation (22) in the non- 
relativistic limit. The function £(y) is an increasing function 
of y. As y — > oo, we have £ — » 1. As y — » a, we have £ — > 0. 
(When 6 = 1, £ is given by eq. A6 in Appendix A.) 

The total energy of the radiation contained in the shock 
front, measured in the frame of the shocked gas, is then 

E R « I (aTi) W^Ar.) « («T 2 4 ) r 3 , (26) 

where the factor 1/3 accounts for the fact that the energy 
density is not uniform (concentrated at the boundary of the 
shock), and the factor 72 in (7 2 Ar s ) accounts for the Lorentz 
contraction. Submitting equation (21) into equation (26), we 
get 



(27) 



47TT S 72 , . 3 2 2 

E R w -2-£,F p (r s v s ) pr F s v s . 

Txl I 

Using the definition of £, we have 

4"/T 7 2C 2^2 2 2„ 

E R ~ -2 F p r T s v s oc r V s v s , 

since F p ~ 1 and 72 /r s ~ 1. 

In the accelerating region (r < -Ra), r s i> s and 72 in- 
crease with r. Hence, the total energy contained in the shock 
front as measured by the rest observer, 72 -Er, increases with 
r. 

In the decelerating region (r > R a , p ~ ?*~ 2 ), by 
equation (13) we have, approximately, r s v s oc r~ 0,2 , thus 
Er oc r 1,8 . In the non-relativistic limit, 72-Bfl, ~ -Bfl oc r 1,8 . 
In the ultra-relativistic limit, 72-E.R oc 72r 1,8 oc r 1,6 since 
7s = r s /\/2 oc r~ 0,2 . Hence, in the region of r > R a , we also 
expect that the total energy contained in the shock front, 
72-Er, increases with r although the shock is decelerating. 
This is caused by the fact that the volume contained in the 
shock front increases with r. 



5 EMERGENCE OF THE SHOCK AND THE 
CHARACTERISTIC QUANTITIES 

Inside the star or deep inside the wind, because of the large 
optical depth in the gas photons have a diffusion velocity 
that is smaller than the velocity of the shock front, so that 
the radiation generated by the shock wave is trapped in- 
side the boundary of the shock. As the shock wave moves 
toward the boundary of the photosphere, the optical depth 
in the gas drops, until a radius is reached where the dif- 
fusion velocity of photons begins to exceed the velocity of 
the shock front. Then, the radiation generated by the shock 
wave starts to escape from the star to produce a bright flash, 
and the shock becomes visible to a remote observer. 

Thus, the shock emerges at a radius where the optical 
depth of the gas to the radiation generated by the shock is 
equal to the optical depth of the shock 

TK = — , (28) 

since beyond that radius photons diffuse outward faster than 
the shock front moves (Imshennik & Nadezhin 1988, 1989; 
Matzner & McKee 1999). Since v s < c always, the shock 
must emerge at a radius where rx > 1. By equations (5) and 
(10), the maximum breakout radius (determined by rx = 1) 
is at 



Vn 



1 - 



1 + — 



1/(1-6) 



(29) 



which is approached by an ultra-relativistic shock. [When 
6=1, the i/max is given by equation (A7).] 

The evolution of T s f3 s is determined by equation (13), 
which can be recast into 



r a /3 s 



v = 



1.181 

xy~° 



E in 



1/2 



1-2 

y 



0.1876 



(30) 



where equations (l)and (2) have been used, and * is the 
mass function defined by equation (16). 

With equations (30), (10), and (5) (or Al if b = 1), 
we can calculate the radius where the shock breakout takes 
place, R br = y br R*, by numerically solving the algebraic 
equation (28). 

After having y br , we can calculate the momentum of the 
shock wave at y = y br , by equation (30). 

Since the shock breakout occurs at a radius where 
t 3 — rx (eq. 28), by equation (27), the total energy of the 
radiation generated by the shock breakout as measured by 
a rest observer is 



[^Er] 



^iF 2 F p pr 3 Y 2 s v 2 a \ 



(31) 



where F p = F p (T s v s ) ~ 1, F 7 = F 7 (r s ) == 7 2 /r s ~ 1. The 
Lorentz factors 72 and F s are related by the equation (5) of 
Blandford & McKee (1976). For the case of 72 = 4/3, F 1 
can be approximated by 



FJx) « — + 



l-l/v 7 ^ 



v/2 ' [1 + 0.9572 (x - 1)]0 9325 > 
which gives the correct asymptotic values at the non- 
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Figure 4. Characteristic quantities of shock emergence as functions of the opacity in the stellar wind. The solid line corresponds to 
e = 10~ 2 . The dashed line corresponds to e = 10~ 3 . Other parameters are: b = 5, E[ n = 10 52 ergs, M e j = IOMq, and i?* = 3Rq. 



relativistic limit (r s — > 1) and the ultra-relativistic limit 
(r s — > oo), and has a fractional error < 0.08% in the trans- 
relativistic case. 

Submitting equations (1) and (2) into equation (31) and 
making use of equation (16), we get 



0.800 x 10 6 K /i ' 25 



„ \ -0.5 

R* \ (ybr\- 

3i? J V 5 ) 



x 1 



-6/4 



(33) 



E hl » *c 2 [^F P T 2 S (3 2 S ] = Vb Jl- — 



1.48 x 10 46 ergs fj, 



y=vb, 
R* 



3Rr. 



i - 



(32) 



The time-duration of the shock breakout event is set 
by the time spent by a photon to diffuse out to the surface 
of the photosphere from the breakout radius. Since at the 
radius of shock breakout the diffusion velocity of photons is 
equal to the velocity of the shock wave, we have 



Similarly, from equation (21), we can obtain the tem- 
perature of the radiation measured in a rest frame 

T br = [ 72 T 2 ] 



( 3 ^° 2 V 4 \ F F i/4 r 3/2.i/2n 



a 



-b/A 



R 



ph,X — Rbr 



R* 



6.96 s 



«s,br A,brC 

R 



(y P h,X - Vbr) 



^ &,br (y P h,X - Vbr) , (34) 



where R p h,x = y P h,xR* is the X-ray photospheric radius 
(eqs. 11 and A4), and v s , hr = Ps,br c = v s( r = -^br) is the 
speed of the shock wave at the time of breakout. 
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Figure 5. Characteristic quantities of shock emergence as functions of the explosion kinetic energy. The solid line corresponds to 
e = 10 -2 . The dashed line corresponds to e = 10~ 3 . Other parameters are: b = 5, k w = 0.7 cm 2 g _1 , M e j = IOMq, and ij* = 3Rq. The 
dotted line shows the solution for shock breakout from a star without a wind, with the same M e j, R+, and re* = 0.2 cm 2 g -1 , £ = 1. 



6 RESULTS 

Unlike in the cases of non-relativistic and ultra-relativistic 
shock waves, where the quantities characterizing the tran- 
sient event from the shock breakout can be expressed 
with factorized scaling relations of input parameters (e.g., 
the eqs. 36-38 of Matzner & McKee 1999), in the trans- 
relativistic case here we must numerically solve the relevant 
equations for the characteristic quantities. 

All the relevant equations are in Sec. 5, supplemented 
by the formulae for the wind mass function ^ and the optical 
depth in Sees. 2 and 3 (and Appendix A when 6 = 1). We 
can eliminate M and «oo from the equations by using 



* = 



807r(b - l)aRi 
k w [(1 - a)i-<> - 1] ' 



(35) 



which is obtained by submitting equation (6) into the defini- 
tion of * (eq. 16). [When b = 1, the corresponding equation 
is (A5).] Since a is a function of e and b (eq. 3), we can then 
choose the input parameters to be Ei n , M c j, i?*, e, b, and 

Note, two opacities are involved in our model: k w , the 
optical opacity in the wind of a WR star, which is used 



to calculate the mass density profile in the wind before the 
shock wave passes through it; kx = k os , the X-ray opacity 
in the wind, which is used to calculate the interaction of the 
X-ray photons generated by the shock wave with particles 
in the wind (see Sec. 2). Since k os = 0.2 cm 2 g _1 is a con- 
stant but somewhat uncertain, we treat k» as an input 
parameter. 

Compared to the case of shock breakout from a star 
without a wind (Matzner & McKee 1999, and Appendix B 
in this paper), here we have two additional parameters: e 
and b, both describing the shape of the wind velocity profile 
(eqs. 2 and 3). However, in the case of a star, the opacity 
is fairly well determined so there are essentially only three 
parameters: the explosion energy E- ln , the ejected mass M e j, 
and the stellar radius i?*. Although there is yet another pa- 
rameter £ = pi/ p* (see Appendix B), which is typically 0.2 
for blue supergiants and 0.5 for red supergiants (Calzavara 
& Matzner 2004), the characteristic quantities (at least the 
energy, the temperature, and the shock momentum) at shock 
breakout are very insensitive to £ (Matzner & McKee 1999). 
While for the problem here, i.e., a dense stellar wind sur- 
rounding a star, the opacity poorly known. Modeling 
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Figure 6. Characteristic quantities of shock emergence as functions of the ejected mass. The solid line corresponds to e = 10 -2 . The 
dashed line corresponds to e = 10~ 3 . Other parameters = 0.7 cm 2 g 1 , _B; n = 10 52 ergs, and ij* = 3Rq. The dotted line 

shows the solution for shock breakout from a star without a wind, with the same E\ n , i?*, and = 0.2 cm 2 g _1 , C = 1- 



of the WR winds indicates that k w is in the range of 0.3-0.9 
cm 2 g _1 at the sonic point, and slightly larger at larger radii 
(Nugis & Lamers 2002). 

The parameter e, which is the ratio of the wind velocity 
at the stellar surface to the terminal velocity of the wind, is 
usually thought to be in the range of 0.001-0.1, and typically 
around 0.01 (Schaerer 1996). 

The parameter b, which characterizes the log-slope of 
the wind velocity in the region near the stellar surface, is 
taken to be unity in the "standard model" of stellar winds. 
However, as already mentioned in Sec. 2, for the winds of 
WR stars 6 can be much larger than unity as argued by 
Robert (1994), Schmutz (1997), and Lepine & Moffat (1999), 
and is typically in the range of 4-6 (Nugis & Lamers 2002). 

In our numerical modeling, we allow k w to vary from 
0.2 to 0.9 cm 2 g -1 , e to vary from 10~ 5 to 10~\ and b from 1 
to 10. We allow Ei n to vary from 10 51 ergs (for normal core- 
collapse supernovae) to 10 53 ergs (for hypernovae), M c j to 
vary from IMq to 2OM0. Although WR stars are compact 
and have small radii, to fully explore the effect of variation 
in the stellar radius on the results, we allow to vary from 
IRq to 3O7?0 . Whenever numbers are quoted, we refer to the 



fiducial values k w = 0.7 cm 2 g \ e = 0.01, 6 = 5, E in = 10 52 
ergs, M e j = 10M Q , and R* = 3Rq, unless otherwise stated. 

Our numerical results for the characteristic quantities of 
the shock breakout, including the total energy (Ebr, eq. 32), 
the temperature (T br , eq. 33), the time-duration (tbr, eq. 34), 
and the shock momentum (r Si br/3s,t>r, eq. 30 with y — j/br) 
are presented in Figs. 4-8. 

Figure 4 shows E hr , T hr , t br , and r Sj brA,br as func- 
tions of the opacity k w . Solid lines correspond to s — 10~ 2 . 
Dashed lines correspond to e = 10~ 3 . Other parameters take 
the fiducial values, as indicated in the figure caption. For 
e = 10~ 2 , E br is a slow but not monotonic function of k w . 
For e = 10~' 5 , Ebr increases with k w . As k w increases from 
0.2 to 0.9 cm 2 g _1 , Sbr increases by a factor ~ 1.2 when 
£ = 10~ 2 , and w 2.6 when e = 10~ 3 . The temperature Tb r 
increases by a factor w 4 in both cases. The opacity k w has 
also an effect on i br , which decreases by a factor when w 4 
when £ — 10~ 2 , and w 2.6 when e = 10 -3 . Similar to the 
temperature, the momentum of the shock front is also an in- 
creasing function of k w . As k w increases from 0.2 to 0.9 cm 2 
g _1 , r 3> biA,br increases by a factor w 2.2 (for both e = 10~ 2 
and £ — 10~ 3 ). Similar to the case of breakout from a star, 



Shock Breakout and GRB 060218/ SN 2006 aj 11 




0.5 1 0.5 1 

log R. (R ) log R. (R ) 

Figure 7. Characteristic quantities of shock emergence as functions of the core radius of the star. The solid line corresponds to e = 10 — 2 . 
The dashed line corresponds to e = 10~ 3 . Other parameters = 0.7 cm g" _1 , E in = 10 52 ergs, and M ej = 1OM . The 

dotted line shows the solution for shock breakout from a star without a wind, with the same Ei n , M e j, and = 0.2 cm 2 g -1 , C = 1- 



the results do not change dramatically with the opacity if b 
is around 5. Thus, the poor knowledge in the opacity in the 
stellar winds will not affect our results drastically. 

All the dependence on k w manifests itself though the 
mass function 9 in equation (35) (and eq. A5 when b = 
1), which shows that oc k^ 1 . Then, from the condition 
for the shock breakout (eq. 28), it can be checked that j/br 
decreases with k w , and (1 — a/j/br)~ & increases with k w . 
From the dependence of £W, Tbr, tbr, and r Sj br/3 s ,br on * 
and j/br, it is not hard to understand the trend shown in 
Fig. 4. First, equation (30) implies that r Sj br/3 s ,br is a strong 
increasing function of k w . Then, equation (34) implies that 
t br is a decreasing function of n w . In equation (32) , * and y br 
decrease with k w , but T 2 bl /3^ br and (1 — a/yhr)~ b increase 
with n w . The overall result on E hr is that shown in the 
top-left panel in Fig. 4. Since the radius of shock breakout 
decreases with k w , the temperature Tb r must increase with 

Figure 5 shows the same set of characteristic quantities 
as functions of the explosion kinetic energy. Symbols and 
values of parameters are similar to Fig. 4 and are explained 
in the figure caption. To compare with the results for a star 



without a wind, we show with dotted lines the corresponding 
characteristic quantities calculated for the shock breakout 
from a star with the formulae in Appendix B, for the same 
values of M c j and R*, and = 0.2 cm 2 g _1 , f = 1. 

As the explosion energy increases from 10 51 ergs to 10 53 
ergs, the breakout energy increases by a factor ~ 117 when 
e = 10~ 2 , and « 188 when e = 10~ 3 . This increasing rate is 
much faster than that in the case of breakout from a star, in 
which the breakout energy increases only by a factor of ~ 22. 
The increase in the temperature is also faster, which is by 
a factor of ~ 33 when s = 10~ 2 , and ~ 55 when e = 10~ 3 
for a star with a dense wind, and only ^ 4.6 for a star 
without a wind. While for the breakout time-duration, it 
appears that for the case of a stellar wind the time-duration 
does not change rapidly when Ei n increases, in contrast to 
the case of a star. This is caused by the fact that when a 
star is surrounded by a dense stellar wind the shock wave 
has more space for acceleration and hence at the time of 
emergence the shock front is more relativistic (see the panel 
for the shock momentum), its velocity approaches the speed 
of light limit. As we have seen in Sec. 5, when the shock 
velocity approaches the speed of light, the breakout radius 
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approaches 2/ m ax- The distance between j/ max and y p h does 
not change with the explosion energy. 

The momentum of the shock front varies with E[ n at 
about the same rate for the case of a stellar wind and the 
case of a star. 

The curvature of the curves in Fig. 5 confirms our claim 
at the beginning of this section that in the trans-relativistic 
case the characteristic quantities of shock breakout in gen- 
eral cannot be written as factorized scaling formulae of input 
parameters. 

Figure 6 shows the characteristic quantities as functions 
of the ejected mass. As the ejected mass M e j increases from 
IMq to 20M Q , the breakout energy decreases by a factor 
w 7.8 when s = 1CP 2 , and ^9.1 when s = 1CP 3 , faster than 
the case of a star for which the decreasing factor £s 4.4. The 
temperature also drops faster. The variation in the breakout 
time-duration is not fast in both the cases of stellar winds 
and stars. That is, the time-duration of the shock breakout 
is not very sensitive to the ejected mass. The momentum of 
the shock front drops slightly slower than that in the case 
of a star. 

Figure 7 shows the characteristic quantities as functions 
of the core radius of the star, which, as in the case of a star 



(Matzner & McKee 1999), is the parameter that most dra- 
matically affects the values of the characteristic quantities. 
As R+ increases from IRq to 3Of?0, the breakout energy 
increases by a factor ~ 69 when e = 1CP 2 , and ~ 51 when 
e = 10 . However, this factor is smaller than that in the 
case of star without a wind, which is w 277. The tempera- 
ture drops very fast, caused by the increase in the area of 
the surface emitting the radiation. As _R* increases from IRq 
to 3Oi?0, the temperature drops by a factor of w 16 when 
e = 10~ 2 , and w 21 when e = 10~ 3 , in contrast to the factor 
w 5.8 in the case of a star. The variation in the stellar radius 
also has a dramatic effect on the breakout time-duration, al- 
though the effect is less prominent than in the case of a star. 
The factor by which the breakout time-duration increases is 
« 43 when e = 1CP 2 , and « 32 when e = 10~ 3 , in contrast 
to that w 590 in the case of a star. The momentum of the 
shock front drops by a factor w 4.8 when e = 10~ 2 , and 
w 4.5 when e = 10~ 3 , similar to the factor w 3 in the case 
of a star. 

Finally, in Fig. 8 we show the characteristic quantities 
as functions of b. The breakout energy E br increases with 
b. As b increases from 1 to 10, E br increases by a factor of 
w 2.9 when e = 10~ 2 , and w 12 when e = 10" 3 . The smaller 
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Figure 9. The observed time-duration of shock breakout (defined by eq. 36) for the models in Figs. 4-8 (from left to right). The solid 
line corresponds to e = 10~ 2 . The dashed line corresponds to e = 10~ 3 . The thin lines are the breakout time-duration tbr without the 
light- travel-time correction (eq. 34). The figure shows that in all models the light-travel-time correction is not important. In some models 
the thin and thick lines are even not visually distinguishable. 



is e, the faster does increases with 6. However, if b is in 
the range of 4-6, the change in Ebr is not essential, only by 
a factor of « 1.2—1.3. The temperature has a similar trend. 
When e = 1CF 2 , the time-duration tbr decreases with b. 
When e = 10 -3 , tbr decreases with b until b grows to a value 
~ 3, beyond which tbr increases but slowly. The momentum 
of the shock front increases with b, caused by the fact that a 
larger b results in a steeper density profile and an enhanced 
acceleration of the shock. 

From Figs. 4-7, and Fig. 8 at b = 4-6, the effects of 
variation in e from 1CF 2 to 1CF 3 can be summarized as fol- 
lows: the breakout energy increases by a factor of 1.8-4; 
the temperature Tb r increases by a factor of 3-5; the shock 
momentum r Sj br/3 s ,br increases by a factor of 2.3-2.5; and 
the time-duration tbr decreases by a factor of 2.3-4.3. 

Figures 5-7 also show that, for a star with a dense wind 
the shock breakout is more energetic than that for a star 
without a wind. This is not surprising, since a star with a 
dense wind has effectively a larger radius so that the shock 
wave has more space and more time for acceleration. For 
the same set of common parameters (E ln , M c - S , R+, but not 
the opacity), for typical parameters the total energy of the 
radiation from shock breakout is larger by a factor > 10 if 
the star is surrounded by a dense wind. The momentum of 
the shock front is also larger by a factor ~ 10. The temper- 
ature does not show a universal trend because of increase in 
the shock breakout radius, but generally it is larger if the 
progenitor is surrounded by a dense wind due to the great 
enhancement in the breakout energy. The time-duration is 
larger for the case of stellar winds as an obvious result of 
increase in the effective radius of the star. 

The shock breakout occurs inside the maximum accel- 
eration radius R a (eq. 14) in all the models presented in 
Figs. 4-8. 

In the above calculations of the breakout time-duration, 
the light-travel-time has not been taken into account. In 
other words, tbr is the duration measured in the supernova 
frame. The duration observed by a remote observer, tbr,obs, 
differs from tbr by an effect caused by the travel-time of 



light — which arises from the fact that an observer will see 
more distant annuli of the stellar disk with a time-delay 
(Ensman & Burrows 1992). The effect of light-travel-time 
could be extremely important when the tbr calculated by 
equation (34) is short, which is definitely true here since 
WR stars are compact. We approximate the observed time- 
duration of the shock breakout event by 



tb 



r,obs 



tt+ti 



ight > 



(36) 



where tn g ht is the light-travel-time. 

In the calculation of the light-travel-time £n g ht, the rela- 
tivistic beaming effect must be taken into account since the 
shock wave in our models is relativistic (Katz 1994). In the 
ultra-relativistic limit, the beaming angle is 9 ~ l/r p h,x, 
where r p h,x is the Lorentz factor of the photosphere which 
can estimated by r ph ,x ~ T s . In the non-relativistic limit, 
we should have 6 = 7r/2. Hence, we use an interpolation 
formula for 9 

(37) 



7r(r s -l) + 2 ■ 
Then, the light-travel-time is 

Rph,X 



flight = 



-(1-COS0) . 



(38) 



When T s > 1, we have t light « R ph ,x/(2F 2 s c). 

For the models presented in Figs. 4-8, we have calcu- 
lated the light-travel-time correction to the observed time- 
duration of the shock breakout event. The results are shown 
in Fig. 9. It turns out that the light-travel-time correction is 
not important. This is caused by the fact that for relativistic 
shock breakout the light-travel-time is significantly reduced 
by the the relativistic beaming effect. 

Numerical results for a set of supernova and WR star 
models are presented in Table 1. From these results we find 
that the efficiency of converting the supernova explosion en- 
ergy to the shock breakout energy, defined by the ratio of 
the breakout energy to the explosion energy, is typically in 
the range of 10 -4 -10 -5 . This efficiency is smaller than that 
in the case of Type II supernovae, which is typically ~ 10~ 3 
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Table 1. Models of Type Ibc supernova explosion and the predicted characteristic parameters for the shock breakout. Input 



parameters: E- ln 


, M e j, R+, 




e, and b. Output parameters: 


?/ph,X, l/bn 


(r s /3 s )br = 




£ br , 


'J- br , ^br.obs 


and fi. 


Models 1- 


4 are 


normal SN 


e Ibc 


Models 


5-7 arc 


hypcrnovac. 














Model 


-C4n 




R* c 


K d 


£ c 


b l 


2/ph,X 


g y br h 






T br k 


^br.obs 




1 


1 


3 


3 


0.7 


0.01 


5 


1.73 


1.45 


1.98 


1.3 


5.4 


2.8 


0.30 


2 


1 


4 


3 


0.2 


0.02 


5 


4.24 


2.45 


0.818 


1.2 


1.9 


25 


1.7 


3 


1.5 


6 


5 


0.5 


0.01 


1 


3.11 


1.61 


0.760 


1.3 


1.7 


35 


2.4 


4 


2 


2 


5 


0.7 


0.001 


5 


1.31 


1.22 


6.72 


19 


28 


1.0 


0.095 


5 


40 


10 


3 


0.2 


0.01 


5 


3.21 


2.53 


5.73 


22 


16 


4.9 


1.0 


6 


50 


10 


10 


0.7 


0.01 


5 


1.73 


1.50 


7.52 


140 


22 


5.5 


0.98 


7 


60 


15 


10 


0.7 


0.002 


5 


1.39 


1.28 


13.7 


320 


62 


2.6 


0.32 



a Explosion kinetic energy in units of 10 51 ergs. 
b Ejected mass in units of Mq. 

°Core radius of the progenitor (the radius at the optical depth of 20), in units of Rq. 
d Optical opacity in the wind, in units of cm 2 g~'. 

° Ratio of the wind velocity at the stellar surface (where r = i?*) to the terminal velocity of the wind (eq. 3). 
' Parameter specifying the profile of the wind velocity (eq. 2). 
s Radius of the X-ray photosphere in units of R±. 

h Radius of the shock front at the time of shock breakout in units of -R*. 
'Momentum of the shock front (eq. 30) at the time of shock breakout. 
JTotal energy of the radiation from the shock breakout, in units of 10 46 ergs. 
k Temperature of the radiation from the shock breakout, in units of 10 6 K = 0.08617 keV. 
'Observed time-duration of the shock breakout event, in units of seconds. 

m Observable defined by the mass-loss rate and the terminal velocity of the stellar wind through eq. (9). 
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Figure 10. Energy versus temperature of shock breakout in Type 
Ibc supernovae produced by core-collapse of a sample of WR stars 
with R pb /R* > 2. 
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Figure 11. Shock momentum versus the observed time-duration 
of shock breakout in Type Ibc supernovae produced by core- 
collapse of the same sample of WRs in Fig. 10. 



if the progenitor is a red supergiant, or ~ 1CP 4 if the pro- 
genitor is a blue supergiant. This is again caused by the fact 
that WR stars have much smaller radii than red and blue 
supergiants. 



7 APPLICATION TO GRB 060218/SN 2006AJ 

As stated in the Introduction, recently it has been claimed 
that supernova shock breakout has been observed in the 
early X-ray emission of GRB 060218, based on the ob- 
servation that a fraction (~ 20%) of the radiation in the 
lightcurve (from 159 s up to ~ 10, 000 s after the trigger 
of the burst) is a soft black-body of temperature w 0.17 
keV (Campana et al. 2006). The total energy estimated for 
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Figure 12. The breakout energy versus the breakout tempera- 
ture (with b = 5). Solid lines correspond to k w = 0.2 cm 2 g _1 . 
Dashed lines correspond to k w = 0.9 cm g . Different solid 
lines (and different dashed lines) correspond to different values 
of e: 10~ 2 , 10~ 3 , 10~ 4 and (upward). Along each line the 

stellar radius varies from IRq (the triangle) to lOOi?0 (the 
point). The supernova explosion energy Ei n = 2 X 10 51 ergs. The 
ejected mass M c j = 2Mq. The region bounded by the dotted line 
indicates the observational constraint on the total energy and the 
temperature of the black-body component in the X-ray afterglow 
of GRB 060218. 

this black-body component is w 10 49 ergs in the 0.3-10 keV 
band, and w 2 x 10 49 in bolometric (S. Campana, private 
communication). A reanalysis carried out by Butler (2006) 
revealed an even larger energy in the black-body, which is 
?s 2 x 10 50 ergs, with a duration w 300 s. 

The overall constraint on the black-body component in 
the early X-ray afterglow of GRB 060218 is summarized as 
follows: the total energy > 10 49 ergs, the temperature is in 
the range of 0.1-0.19 keV (i.e., 1.2-2.2 x 10 6 K), and the 
duration > 300 s (Campana et al. 2006; Butler 2006). 

In this Section, we apply the procedure developed in 
previous sections to calculate the characteristic quantities of 
the shock breakout event for SN 2006aj with the assumption 
that the supernova was produced by the core-collapse of a 
WR star surrounded by a dense wind, and examine if the 
black-body component in GRB 060218 can be interpreted 
by the shock breakout in SN 2006aj. 

First, we apply the procedure to the WR stars in Fig. 1, 
which are among the best studied catalog of WR stars with 
model-determined stellar and photospheric radii (Hamann 
et al. 1995; Koesterke & Hamann 1995; Grafener et al. 1998). 
We pick up only stars with j/ p h = R p h/R* > 2, since oth- 
erwise the e given by our simplified model would be too 
small (Fig. 3). The sample so selected consists of total 36 
stars, including 20 Galactic WCs, 10 Galactic WNs, and 6 
LMC WCs. The majority of WC stars have been included. 
Since theses stars were modeled by the standard stellar wind 
model (Appendix A), we choose 6=1. The a parameter is 
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Figure 13. Same as Fig. 12 but with 6=1. 



then obtained from the published values of y p h = R P h/R* 
in the papers cited above by equations (A2) and (A3). The 
value of e is calculated with equation (3). Then, with the 
published values of M, Voo , and i?*, the constant mean opac- 
ity k w can be calculated with equation (A5). 

For the supernova explosion energy and the ejected 
mass, we take the values obtained by modeling the spec- 
tra and the lightcurve of SN 2006aj: E in = 2 x 10 51 ergs, 
and M cj = 2M (Mazzali et al. 2006). Then, for each star 
we have all of the six parameters needed for calculating the 
characteristic quantities of shock breakout. 

Our results are shown in Fig. 10 for the breakout en- 
ergy versus the breakout temperature, and in Fig. 11 for 
the breakout shock momentum versus the breakout time- 
duration. Note, here the breakout time-duration has in- 
cluded the light-travel-time (eq. 36), so it corresponds to 
the observed time-duration. From Fig. 10 we see that, al- 
though the temperature is in the range of the black-body 
component in GRB 060218 for several WC stars (on the left 
end) , the total energy of the radiation arising from the shock 
breakout never exceeds 10 47 ergs, i.e., always smaller than 
the total energy of the observed black-body component in 
GRB 060218 by more than two orders of magnitude. 

From Fig. 11, the time-duration of the shock breakout 
never exceeds 100 s, also well below the observational limit 
on the black-body component in GRB 060218. 

Thus, it appears that none of the stars in the considered 
sample of WRs is able to produce a supernova with shock 
breakout energy that is large enough to explain the black- 
body component observed in the early X-ray afterglow of 
GRB 060218. 

Of course, GRBs are rare events compared to super- 
novae (Podsiadlowski et al. 2004) , they may require progen- 
itors that are in more extreme conditions than the WRs 
in our sample. To test the possibility for explaining the 
black-body component in GRB 060218 with shock break- 
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out from WR stars in a larger parameter space, we have 
calculated a number of models with a large range of param- 
eters. The results are shown in Fig. 12 (b — 5) and Fig. 13 
(b — 1). The explosion energy and the ejected mass are fixed 
at Ei n = 2 x 10 51 ergs and M c j = 2M@ , as obtained by mod- 
eling the spectra and the lightcurve of SN 2006aj (Mazzali 
et al. 2006). We allow e to vary from 10~ 2 to 10~ 5 . For 
the opacity k w , we choose two extreme values: 0.2 cm 2 g _1 
(solid lines) and 0.9 cm 2 g _1 (dashed lines). The observa- 
tional bound on the total energy and the temperature of the 
black-body component in the early X-ray emission of GRB 
060218 is shown in the figures by the region bounded by the 
dotted lines. 

The radius of the star, i?*, which is the parameter that 
the characteristic quantities of the shock breakout are most 
sensitive to, is allowed to vary from IRq to 100/?©, covering 
a space of radii that is more than enough for WR stars. 

Figures 12 and 13 show that to explain the black- 
body component observed in the early X-ray emission of 
GRB 060218, the radius of the progenitor WR star must 
be > lOOi?0. It is very unlikely that there exist WR stars 
having so large stellar radii. Although it is possible to get 
Ebr > 10 49 ergs with i?* < 100i?o if e is very small and/or 
k w is very large, the corresponding Tb r would be too high to 
be consistent with the temperature of the black-body com- 
ponent in GRB 060218. 



8 SUMMARY, CONCLUSIONS, AND 
DISCUSSIONS 

We have presented a simple model for calculating the charac- 
teristic quantities (total energy, temperature, time-duration, 
and shock momentum) for the flashes arising from shock 
breakout in Type Ibc supernovae produced by the core- 
collapse of Wolf-Rayet stars surrounded by dense stellar 
winds. The wind velocity is modeled by equation (2) , a pro- 
file that is often adopted in the study of stellar winds. How- 
ever, in contrast to the case for O-stars where the parameter 
6 is close to unity, for WR star winds 6 can be much larger 
and is usually in the range of 4-6 (Nugis & Lamers 2002). 
The opacity in the wind, assumed to be a constant, 

which is a reasonable approximation for the calculation of 
the optical depth since the opacity varies with radius very 
slowly compared to the mass density of the wind (Nugis & 
Lamers 2002). Modeling of the opacity in the winds of WR 
stars indicates that k w is in the range of 0.3-0.9 cm 2 g _1 
(Nugis & Lamers 2002). 

Our model is an extension of the existing model for 
calculating the characteristic quantities for supernova shock 
breakout from a star without a wind, which is suitable for 
Type II supernovae (Imshennik & Nadezhin 1988, 1989; 
Matzner & McKee 1999; Tan et al. 2001). Due to the com- 
pactness of WR stars, the shock momentum is expected to 
be trans-relativistic at the time of breakout. Thus, we have 
followed Blandford & McKee (1976) and Tan et al. (2001) 
to take into account the relativistic effects. 

Because of the large optical depth in the wind, the su- 
pernova shock breakout occurs in the wind region rather 
than in the interior of the star. This is equivalent to say that 
the presence of a dense stellar wind effectively increases the 
radius of the star. As a result, the shock has more space and 



more time for acceleration, and the shock breakout appears 
to be more energetic than in the case for the same star but 
the effect of the stellar winds is not taken into account (see. 
e.g., Blinnikov et al. 2002). 

The formulae for determining the radius where the 
shock breakout occurs and that for computing the charac- 
teristic quantities for the radiation arising from the shock 
breakout are collected in Sec. 5. They include equations (28), 
determining the breakout radius; (30), evaluating the mo- 
mentum of the shock; (32), (33), and (34), calculating the 
energy, temperature, and the time-duration of the radiation 
from shock breakout. Although exact and analytic solutions 
are impossible because of the trans-relativistic nature of the 
problem, all the equations are algebraic and a simple nu- 
merical program is able to calculate all the characteristic 
quantities. The model contains six input parameters: the 
explosion kinetic energy (£"i n ), the ejected mass (M e j), the 
core radius of the star (R+, the radius where the optical 
depth t w — 20), the opacity in the wind (k w ), the param- 
eter b specifying the wind velocity profile, and the ratio of 
the wind velocity at the stellar surface (where r = R+) to 
the terminal velocity of the wind (e). 

Our numerical results are summarized in Figs. 4-8 and 
Table 1. Figs. 4-8 illustrate how the characteristic quantities 
vary with the input parameters. As in the case of shock 
breakout from a star without a wind, the core radius of the 
star is the most important parameter affecting the results. 
That is, the characteristic quantities are most sensitive to 
the variation in the stellar radius. This feature leads to the 
possibility for distinguishing the progenitors of supernovae 
by observing the flashes from the shock breakout (Calzavara 
& Matzner 2004). In addition, in the case of dense stellar 
winds, the results are more sensitive to the variation in the 
supernova explosion kinetic energy. For example, roughly 
speaking, Ebr oc E ln when the star has a dense wind, in 
contrast to Ebr oc -E ; ° n 6 in the case of a star without a wind. 
Overall, the shock breakout from a star with a dense wind 
is more energetic than that from a star without a wind. 
For a star of the same radius, and for the same explosion 
kinetic energy and ejected mass, the total energy released 
by the shock breakout is larger by a factor > 10 if the star 
is surrounded by a thick wind. The time-duration is also 
larger, and the shock momentum at the time of breakout is 
more relativistic. 

For explosion energy Ei n — 10 51 ergs, ejected mass 
M e j = 3Mq , and stellar radius _R* = 3Rq (typical values for 
normal SNe Ibc), we get breakout energy Ebr ~ 1.3 x 10 46 
ergs, temperature T br « 5.4 x 10 6 K « 0.46 keV, and ob- 
served time-duration tbr.obs ~ 2.8 s if other parameters take 
fiducial values (k = 0.7 cm g~\ b = 5, and e = 0.01). For 
E in = 5 x 10 52 ergs, M oj = 10M©, and ii* = 10i?© (typ- 
ical values for hypernovae), we get Ebr w 1.4 x 10 48 ergs, 
Tbr « 2.2 x 10 7 K « 1.9 keV, and t br ,obs « 5.5 s. More 
numerical results are shown in Table 1. 

We have applied our model to GRB 060218/SN 2006aj, 
in which a soft black-body component has been observed 
in the early X-ray emission of the GRB and has been in- 
terpreted as an evidence for the supernova shock breakout 
(Campana et al. 2006). We take the values of the supernova 
explosion energy and the ejected mass obtained by modeling 
the spectra and the lightcurve of the supernova (Mazzali et 
al. 2006). We find that, the energy released by the super- 
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nova shock breakout in a thick wind of a WR progenitor 
star is generally too small to explain the black-body radi- 
ation in GRB 060218. To obtain the breakout energy and 
the temperature that are consistent with the observational 
constraint, the core radius of the progenitor WR star has 
to be > 100i?©, which is much too large for a WR star. 
Thus, we conclude that the black-body component in the 
X-ray afterglow of GRB 060218 cannot be interpreted by 
the shock breakout in the underlying supernova. Instead, it 
must originate from other processes which might be related 
to the GRB outflow (see, e.g., Fan, Piran & Xu 2006). This 
conclusion is in agreement with the analysis by Ghisellini, 
Ghirlanda & Tavecchio (2006). 

One may argue that GRB-connected supernovae should 
be highly aspherical so that our spherical model might have 
under-estimated the energy of the shock breakout. The ef- 
fect of explosion asymmetry can be estimated as follow. 
Assume that the explosion produces a shock wave in a 
solid angle ft = 4-kuj < 4% with a kinetic explosion en- 
ergy Ei n , which ejects a mass M e j from the progenitor. The 
shock wave is symmetric in the azimuthal direction and 
does not expand to the outside of fi. The motion of the 
shock wave would then be the same as that of a spher- 
ical shock wave (ui = 1) with a kinetic explosion energy 
uj~ 1 Ei n and an ejected mass u)~ 1 M e j, assuming that the 
progenitor is spherically symmetric. Then, by equation (30), 

p oc < /2 M c 7 - 313 ^- - 187 = (^-°- 374 i5 in ) 1/2 M c 7 - 313 . That 
is, the motion of the asymmetric shock wave can be calcu- 
lated by equation (30) but with E[ n replaced by a larger 
E' in = cu' - 374 E in . Then, by Fig. 5, the temperature T br , the 
shock momentum F Sj br/3s,br, and the isotropic- equivalent en- 
ergy Ebr of the asymmetric shock breakout are larger than 
that in a spherical explosion with the same E- 1Yl and M e j. 
However, the time-duration t br is not sensitive to lo. 

Indeed, aspherical explosion has been claimed to be ob- 
served in the luminous Type Ic SN 2003jd, in which the 
double-lined profiles in the nebular lines of neutral oxygen 
and magnesium revealed in later-time observations by Sub- 
aru and Keck are explained as results of observing an as- 
pherical supernovae along a direction almost perpendicular 
to the axis of the explosion (Mazzali et al. 2005). However, 
for SN 2006aj, there is no any evidence for aspherical ex- 
plosion. Observation on the radio afterglow and modeling of 
it indicate that the outflow associated with GRB 060218 is 
mildly relativistic so should be more or less spherical (Soder- 
berg et al. 2006; Fan, Piran & Xu 2006, see also Li 2006). 

We should also remark that whether the progenitors 
of GRBs are surrounded by dense winds is still an open 
question. Although a wind-type density profile is naturally 
expected for the environment surrounding a GRB as its pro- 
genitor is broadly thought to be a massive star, observations 
on the GRB afterglows have revealed that most of the af- 
terglow data are consistent with a constant density external 
medium and only a handful of bursts can be well modeled by 
the wind model (Berger et al. 2003; Zhang & Meszaros 2004; 
Panaitescu 2005; Fryer, Rockefeller & Young 2006, and ref- 
erences therein). For the case of GRB 060218, modeling of 
its radio afterglow also does not favor a dense circum-burst 
wind profile (Soderberg et al. 2006; Fan et al. 2006). 

A theoretical argument against strong winds surround- 
ing GRB progenitors comes from the consideration of an- 



gular momentum (Yoon & Langer 2005; Woosley & Heger 
2006a, and references therein) . For a black hole formed from 
the core-collapse of a massive star to have a disk rotating 
around it and to launch a relativistic jet, the progenitor star 
must rotate rapidly with the specific angular momentum in 
the core j > 3 x 10 16 cm 2 s" 1 (MacFadyen & Woosley 1999). 
To satisfy this requirement, the progenitor star should not 
have had a phase with an intense stellar wind since a dense 
wind is very effective in removing angular momentum. Given 
the fact that the mass-loss rate of a star sensitively depends 
on its metallicity (Vink & de Koter 2005) and the observa- 
tions that GRBs prefer to occur in galaxies with low metal- 
licity (Fynbo et al. 2003; Hjorth et al. 2003a; Le Floc'h et 
al. 2003; Sollerman et al. 2005; Fruchter et al. 2006; Stanek 
et al. 2006), it is reasonable to expect that the progenitors 
of GRBs should not have dense stellar winds surrounding 
them. Even in this situation, however, the radius of the mas- 
sive progenitor star is also very unlikely to be large enough 
(> IOO-R0) to explain the black-body component in GRB 
060218 since its progenitor star has only a mass ~ 20M© as 
obtained by modeling the supernova lightcurve and spectra 
(Mazzali et al. 2006). In addition, if the progenitor does not 
have a thick wind, then in calculating the results for the 
shock breakout one should use the formulae in Appendix B 
for a star without a wind. But in Sec. 6 we have seen that 
the formulae for a star without a wind lead to smaller total 
energy in the radiation from the shock breakout than the 
formulae for a star with a dense wind. 

In spite of the disappointing result on GRB 060218/SN 
2006aj, our model is expected to have important applica- 
tions to Type Ibc supernovae since whose progenitors arc 
broadly believed to be WR stars. In addition, some Type 
II supernovae appear also to be related to progenitor stars 
with intensive stellar winds, e.g. SNc Iln (also called Ildw) 
(Hamuy 2004). Observations on the transient events from 
supernova shock breakout will be the most powerful ap- 
proach for diagnosing the progenitors of supernovae. For 
this goal we would like to mention LOBSTER, an upcom- 
ing space observatory dedicated to detect soft X-ray flashes 
from shock breakout in supernovae (Calzavara & Matzner 
2004). 
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APPENDIX A: OPTICAL DEPTH IN A WIND 
IN THE STANDARD MODEL 

In the standard model of stellar winds the parameter b in 
equation (2) is assumed to be unity. Then, the integral in 
equation (4) gives 



t w = T In 1 



where 



20 



(Al) 



(A2) 



aR* ln(l-a) -1 ' 

The ratio of the photospheric radius (at t w = 2/3) to 
the stellar core radius (at t w = 20) is 
a 



l-exp[-2/(3r )] ' 

which approaches 1 as a — » 1, and 30 as a — > 0. 
The X-ray photospheric radius is at 



y P h,x 



1 - exp[-2/(3tro)] 



(A3) 



(A4) 



The corresponding mass function * (eq. 16), when M 
and Voo are eliminated (by using eq. A2) , is 



SOnaRi 



k w ln(l — a) 1 

The parameter £ = IdhiTx/dlnrp 1 (Sec. 4) is 



a 



V 



(A5) 



(A6) 



which approaches unity as y — > oo, and approaches zero as 
V -> a. 

The maximum radius where the shock breakout occurs 
(see Sec. 5) is given by 
a 



1 — exp(— 1/tro) 



(A7) 



APPENDIX B: SHOCK BREAKOUT FROM A 
STAR WITHOUT A WIND 

The mass density in an outer layer of a star is described by 
a power law (see, e.g., Matzner & McKee 1999) 



P = Pix , 



(Bl) 



where x = 1 — r/R+, n is related to the polytropic index 7 
by 7 = 1 + 1/n. When 7 = 4/3, we have n = 3. 
The optical depth in the star is 

n+l _ _ n*piR* 



T X 



To 



(B2) 



n + l ' 

where k* is the opacity. 

Near the stellar surface we have r w i?*, so the shock 
accelerates according to equation (13) with m w M c j and 
r w ii* (Tan et al. 2001). 
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Figure CI. The mass function defined by eq. (16) against the 
photospheric radius for the sample of WRs in Fig. 1. Clearly there 
is a strong correlation between ^ and Rph- The solid straight line 
is the best fit to all the data by eq. (CI). The dashed straight line 
is the best fit to the WC stars (filled symbols) by eq. (C2). 



The geometric thickness of the shock front is 



Ar s 



r s R+x 



(B3) 



r 2 sK *p V* ri ' 

where £ = l/(n + 1), t s = c/v s . 

The shock breakout occurs at a radius where r* = t s . 
The minimum value of x^, which occurs when v s — > c, is 

x min = r- 1/(n+1) , (B4) 

corresponding to the maximum breakout radius r max = 

-R*(l — Xmin). 

The pressure of the gas behind the shock front, mea- 
sured in the frame of the shocked gas, is still given by equa- 
tion (21), from which the temperature of the shock emer- 
gence can be calculated. 

The total energy of radiation in the shock emergence, 
measured in the rest frame, is 



E hl w 4TY^F p pRt(r s v s fx\ r ^ R ^ . 

The time-duration of the shock breakout is 
R*x hl 



tbr 



(B5) 



(B6) 



The input parameters include _Ei n , M e j, R±, k+, and 
£ = pi/p*, where p* = M e JR+. 



APPENDIX C: A CORRELATION IN 
WOLF-RAYET STAR PARAMETERS 

From the parameters of the 92 Galactic and LMC WR stars 
presented in Fig. 1, a correlation between \& = MR* /voo 
(eq. 16) and 7? ph can be derived. 

In Fig. CI, we plot log* against logi? p h for the 92 



20 Li-Xin Li 



WRs. Clearly, there is a strong correlation between * and 
i?ph. The relation is best fitted by 



log* = -10.87 + 2.36 logi? ph 

for all stars, and 

log* = -10.42 + 1.77 logi? ph 



(CI) 
(C2) 



for WC stars only, where * is in units of Mq , and R p h is in 
units of Rq. 

To the knowledge of the author the relation does not 
exist in the literature so is presented here, although it is 
irrelevant to the subject of the paper. 

This paper has been typeset from a TjnjX/ WJ^K file prepared 
by the author. 



